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Abstract. We present pointwise gradient bounds for solutions to p-Laplacean 
type non-homogeneous equations employing non-linear Wolff type potentials, 
and then prove similar bounds, via suitable caloric potentials, for solutions to 
parabolic equations. A method of proof entails a family of non-local Cacciop- 
poli inequalities, together with a DeGiorgi's type fractional iteration. 



1. The CLASSICAL SETTING AND A ZERO ORDER ESTIMATE 

In this note we describe some of the results and techniques developed in the 
papers O [12] , which give a complete non-linear analog of the classical pointwise 
gradient estimates valid for the Poisson equation 

(1.1) -Au^fi inM", 

where /i is in the most general case a Radon measure with finite total mass. More- 
over, the estimates we present hold for non-linear parabolic equations. At the same 
time our results give a somehow unexpected but natural maximal order - and par- 
abolic - version of a by now classical result due to Kilpelainen & Maly [T7] and 
later extended, by mean of a different approach, by Trudinger & Wang [24j . To 
better frame our setting, let us recall a few basic linear results concerning the basic 
example (|l.ip - here for simplicity considered in the whole M" - for which, due to 
the use of classical representation formulas, it is possible to get pointwise bounds 
for solutions via the use of Riczs potentials 

(1.2) ^^(0'"] 
such as 

(1.3) \u{x)\ < chilfiDix) , and \Du{x)\ < cIli\^i\)ix) . 

We recall that the equivalent, localized version of the Riesz potential Ip{ii){x) is 
given by the linear potential 

(1.4) lp[xo,R) :== / , /? e (0,n] 

Jo Q ' Q 

with B(xo, q) being the open ball centered at xq, with radius g. In fact, it is not 
difficult to see that 

(1.5) I';ixo,R)< f -^^^^^^Ip{p.^B{xo,R)){xo)<Up){xo) 

Jbr{x„) fo -y\ 

holds provided /i is a non-negative measure. A question is now, is it possible to give 
an analogue of estimates ()1.3p in the case of general quasilinear equations such as 
for instance, the degenerate p-Laplacean equation 



(1.6) ^ div {\Du\P-^ Du) ^ 



? 
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A first answer has been given in the papers [T71 [H] , where - for suitably defined 
solutions to ()1.6p - the authors prove the following pointwise zero order estimate - 
i.e. for u - when p < n, via non-linear Wolff potentials: 

(1.7) \uixo)\<c{~f \u\P-Ux] +cW'^,^{xo,2R), 

\Jb{xo,R) ) 

where the constant c depends on the quantities n, p, and 

(1.8) W^p(.To,i?) :=y^ j — /3e(0,7i/p], 

is the non-linear Wolff potential of /i. Of course we are here using the standard 
notation concerning integral averages 

/ \u\'^ dx -.^ ^ / \u\'^ dx . 

\B{xo,R)\ 

JB{xa,R) 

Estimate (|1.7p . which extends to a whole family of general quasi-linear equations, 
and which is commonly considered as a basic result in the theory of quasi-linear 
equations, is the natural non-linear analogue of the first linear estimate appearing 
in p.3p . Here we present the non-linear analogue of the second estimate in 
thereby giving a pointwise gradient estimate via non-linear potentials which upgrades 
(|1.8p up to the gradient /maximal level. 



2. Degenerate Elliptic estimates 

In this section the growth exponent p will be a number such that p > 2, we shall 
therefore treat possibly degenerate elliptic equations when p 2. Specifically, we 
shall consider general non-linear, possibly degenerate equations with p-growth of 
the type 

(2.1) — div a(.T, Du) = /i . 

whenever /Lt is a Radon measure with finite total mass defined on eventually 
letting /x(R" \ J7) = 0, without loss of generality we may assume that ^ is defined 
on the whole M". The continuous vector field a: f2 x R" R" is assumed to be 
C^-regular in the gradient variable z, with az{-) being Carathcodory regular and 
satisfying the following growth, ellipticity and continuity assumptions: 

{ |a(x,z)| + |a,(x,z)|(|z|2 + s2)^ <L(|z|2 + s2)^ 

(2.2) J i^-'d^P + s2)^ |A|2 < {a,{x, z)A, A) 

[ \a{x, z) — a{xQ, z)\ < LiU!{\x - xa\){\z\'^ + s^)^ , 

whenever x,xq G Vt and z,A G R", where < < 1 < L and s > 0,-Li > 1 
are fixed parameters. When p > 2 we shall assume that there exists a positive 
a < min{l,p — 2} such that the Holder continuity property 

(2.3) \a,{x, Z2) - a,(x,zi)| < L\z2 ~ ^iTd^iP + 1^21' + s^)^^ 

holds whenever zi,Z2 G R" and x G ^l. Here uj: [0, oo) — > [0, cxd) is a modulus of 
continuity i.e. a non-decreasing function such that u)(0) = and lu{-) < 1. On such 
a function we impose a natural decay property, which is essentially optimal for the 
result we are going to have, and prescribes a Dini continuous dependence of the 
partial map x i— > a{x, z)/(\z\ -\- s)p~^: 

(2.4) rHg)]-p^:^d{R)<^, 

Jo Q 
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for some R > 0. The prototype of (|2.ip is - choosing s = and omitting the x- 
dependence - clearly given by the p-Laplacean equation ()1.6p . which satisfies ()2.3p 
whenever a < min{l,j3 — 2}. In the following, when a measure /x actually turns out 
to be an L^-function, wc shall use the standard notation 



whenever A is a measurable set on which fi is defined. 

In this paper we shall present our results in the form of a priori estimates - i.e. 
when solutions and data are taken to be more regular than needed, for instance 
u 6 C^(r2) and fi G L^(il) - but they actually hold, via a standard approximation 
argument, for general weak and very weak solutions - i.e. distributional solutions 
which are not in the natural space W^^'^(57) - to measure data problems such as, 
for instance 



where /i is a general Radon measure with finite total mass, defined on Q. The 
reason for such a choice is that the approximation argument in question leads to 
different notions of solutions, according to the rcgularity/integrability properties of 
the right hand side fi. We do not want to enter in such details too much, for which 
we refer to }12[ I22j , and therefore we confine ourselves to the neat a priori estimate 
form of the results. 

For instance, in the case (j2.5p with fi being a genuinely Radon measure, in 
]TI[ [22] we consider the so called Solutions Obtained by Limit of Approximations 
(SOLA), which is a standard class considered when dealing with measure data 
problems. Such solutions are in particular unique in the case p = 2, as proved in 
[6l [25] . Finally, if the right hand side of (|2.1|) is integrable enough to deduce that 
fi e W~^'P (f2), then our results apply to general weak energy solutions u G W^'^{il.) 
to (fTT|) . 

The first result we present is now 

Theorem 2.1 (Non-linear potential gradient bound). Let u £ C^(f2), he a weak 
solution to (j2.ip with ^ G L^{il), under the assumptions (j2.2p . Then there exists 
a constant c = c{n,p,i/, L, Li,a) > 1, and a positive radius Rq depending only on 
n,p,v^ L, Li^Lj{-),a, such that the pointwise estimate 



holds whenever B{xo,2R) C fl^ and R < Rq. Moreover, when the vector field a(-) 
is independent of x - and in particular for the p-Laplacean operator p. 611 - estimate 
(|2.6p holds with no restriction on R. 

The potential appearing in (|2.6|) is the natural one since its shape respects 

the scaling properties of the equation with respect to the estimate in question; 
compare with the linear estimates (|1.3p . When extended to general weak solutions 
estimate (|2.6p tells us the remarkable fact that the boundedness of Du at a point 
Xq is independent of the solution u, and of the vector field a(-) considered, but only 
depends on the behavior of in a neighborhood of xq. 

A particularly interesting situation occurs in the case p = 2, when we have a 
pointwise potential estimate which is completely similar to the second one in (jl.3p . 
and that we think deserves a statement of its own, that is 




(2.5) 




in Q 
on dil 



(2.6) 
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Theorem 2.2 (Linear potential gradient bound). Let u £ C^{^1), be a weak 
solution to ()2.ip with ^ g L^(n), under the assumptions (|2.2p considered with 
p = 2. Then there exists a constant c = c{n,p,h', L, Li) > 0, and a positive radius 
Rq = Ro{n,p,i', L, Li,Lu(-)) such that the pointwise estimate 



(2.7) 



\Du(xq)\ < c 



/ 



R) 



{\Du\ 



I dx 



holds whenever B{xq,2R) C $7, and R < Rq- Moreover, when the vector field a(-) 
is independent of the variable x, estimate (|2.7p holds with no restriction on R. 



Beside their intrinsic theoretical interest, the point in estimates ()2.6p - (|2.7p is 
that they allow to unify and recast essentially all the gradient L'-estimates for 
quasilinear equations in divergence form; moreover they allow for an immediate 
derivation of estimates in intermediate spaces such as interpolation spaces. We 
refer to the recent survey [Hj for an account of such estimates. Indeed, by (|2.6p it 
is clear that the behavior of Du can be controlled by that , which is in turn 

known via the behavior of Riesz potentials. In fact, this is a consequence of the 
pointwise bound of the Wolff potential via the Havin-Maz'ja non linear potential 
[HIHIEI, that is 



(2. 



(-,00) 



— < cli 

Q P 



(a;o)- 



Ultimately, thanks to 
we have 



and to the well-known properties of the Riesz potentials. 



(2.9) 



n<l(p-l) 



q e (l,n), 



while Marcikiewicz spaces must be introduced for the borderline case q = 1. In- 
equality (|2.9p immediately allows to recast the classical gradient estimates for so- 
lutions to (|2.5p such as those due to Boccardo & Gallouet [II |8] - when q is "small" 
- and Iwaniec [16] and DiBenedetto & Manfredi [10] - when q is "large" - that is, 
for solutions to (12.51) it holds that 



^ e L'' =^ Du e L "-9 ge(l,rt). 

Moreover, since the operator fi t-^ is obviously sub-linear, using the estimates 

related to (|2.9p and classical interpolation theorems for sub-linear operators one 
immediately gets estimates in refined scales of spaces such Lorentz or Orlicz spaces, 
recovering some estimates of Talenti [23], but directly for the gradient of solutions, 
rather than for solutions themselves. 

Another point of Theorem 12.11 is that it allows to prove an essentially optimal 
Lipschitz continuity critcrium with respect to the regularity of coefficients (|2.4p . 
that is 



(2.10) 



(•, R) e L°°(f7), for some R > =^ Du e L^,{n, 



and moreover the local bound 



(2.11) \\Du\\ 



holds whenever B2R C f2 



< c 



'B{xo,R) 



\Du\+s)^dx I +c Wt^{-,R) 
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We finally recall that another consequence of the classical estimate 



and of 



IS 



(2.12) \Du{xo)\<c{f {\Du\ + s)idx] +cli\ /i(|/i|) 



which holds whenever B[xo,2R) C satisfies the conditions imposed in Theorem 
12.11 Here we recall the reader that we have previously extended fu, to the whole 
space M". 



3. Parabolic first, and zero order estimates 

Our aim here is not only to give a parabolic version of the elliptic estimate ()2.6p , 
but also to give a zero order estimate, that is the parabolic analog of the zero order 
elliptic estimate jl7] . the validity of which was yet considered to be an open issue. 
We consider quasilinear parabolic equations of the type 

(3.1) ut — div a(.T, t, Du) ~ /i , 

in a cylindrical domain CIt := x (— T, 0), where as in the previous section Q C 
R", n > 2 and T > 0. The vector-field a: f^T x R" R" is assumed to be 
Caratheodory regular together with az{-), and indeed being C^-rcgular with respect 
to the gradient variable z G R", and satisfying the following standard growth, 
ellipticity/parabolicity and continuity conditions: 



(3.2) 



\a{x,t,z)\ + \a^{x,t,z)\{\z\ + s) < L{\z\ + s) 
< {a,{x,t,z)\\) 
\a[x, t, z) — a{xo,t, z)\ < LiU}{\x - a;o|)(|z| + s 



for every choice of x, xq G fi, z, A G R" and t G {—T, 0); here the function oj : R+ ~> 
R"*" is as in (PUJa. Note that anyway we are assuming no continuity on the map 
1 1— > a{-,t, •), which is considered to be a priori only measurable. In other words we 
are considering the analog of assumptions (|2.2p for p = 2; the reason we are adopting 
this restriction is that when dealing with the evolutionary p-Laplacean operator 
estimates assume the usual form only when using so called "intrinsic cylinders", 
according the parabolic p-Laplacean theory developed by DiBenedetto [H] . These 
are - unless p = 2 when they reduce to the standard parabolic ones - cylinders whose 
size locally depends on the size of the solutions itself, therefore a formulation of 
the estimates via non-linear potentials - whose definition is built essentially using a 
standard family of balls and it is therefore "universal" - is not immediate and will 
be the object of future investigation. We refer to [I] for global gradient estimates. 

In order to state our results we need some additional terminology. Let us recall 
that given points {x,t), {xo,to) G R""*"^ the standard parabolic metric is defined by 



(3.3) dpar((a;,t), (a;o,to)) := max{|a; - xq\, V|t - toll « Vl^ ~ ^o|^ + 1^ - *o| 

and the related metric balls with radius R with respect to this metric are given by 
cylinders B{xo,R) x {to — R'^,to + R^). The "caloric" Ricsz potential - compare 
with elliptic one defined in (|1.2p . and with [2], for instance - is now built starting 
from (ESI) 



(3.4) /,(A^)((x,t)):= / rrlw'lU-/^ ' 0<(3<N:^n + 2, 

JR" + i dpar((x,i), (X^tjy P 

whenever (x,t) G R"+^. In order to be used in estimates for parabolic equations, 
it is convenient to introduce its local version via the usual backward parabolic 
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cylinders - with "vertex" at {xo,to) - that is 

(3.5) Q{xo,to;R) ■.^B{xo,R) x {to ^ R^to) , 
so that wc define 

(3.6) i;;(xo, to; R) £ ^^iQM■, 9)) | ^^^^^ ^ ^ N] . 
The main result in the paraboUc case is 

Theorem 3.1 (ParaboHc potential gradient bound). Under the assumptions p.2p 
and (1231), let u G C°(-T, 0; ^^(f^)) &e a weak solution to with fi G L°°i^T) 

and such that Du G C'^(IIt)- Then there exists a constant c = c{n,v,L) and a 
radius Rq = RQiji^v, L, Li,lu{-)) such that the following estimate: 



■JQ{xo,ta; 



(3.7) \Du{xo,to)\ < c~f {\Du\+s)dxdt + cl[^^{xo,to;2R). 



holds whenever Q{xQ,tQ;2R) C fl, and R < Rq. When the vector field a(-) is 
independent of the space variable x, estimate ()3.7p holds with no restriction on R. 



Again, as in the eUiptic case, estimate (|3.7p also holds for solutions to general 
measure data problems as 



(3.8) 



ut — diva(a;, t, Du) = fi in ^It 
u = on dpai-flr 



where /x is a general Radon measure with finite mass on ^It, that we shall again 
consider to be defined in the whole R"+^. In the spirit of the elliptic result (j2.1ip 
we have the following implication, which provides a boundedness criteria for the 
spatial gradient, under the Dini continuity assumption for the spatial coefficients 
stated in ([^ : 

(3.9) Ii^'(-;^) e L°°{nT), for some R > =^ Du e LfJe(f7T,K") . 

We conclude with the zero order potential estimate, which applies to general equa- 
tions of the type (|3.ip when considered with a measurable dependence upon the 
coefficients {x,t). The relevant hypotheses here are the following standard growth 
and monotonicity properties: 



(3.10) 



|a(x, t, z)\ < L{\z\ + s) 
v\z2 - ^iP < {a{x, t, Z2) - a{x,t, zi),Z2 - zi) 



which are assumed to hold whenever {x,t) G fir and z,zi,Z2 G M". In particular, 
since the pointwise bound will be derived on m, rather than on Du, wc do not need 
any differentiability assumption on a(-) with respect to the spatial gradient variable 
z-variable, assumptions p.lOp are clearly weaker than 



Theorem 3.2. Under the assumptions dSHU]), letu G L^{-T,0;W^-^{n))r\C° (nr) 
be a weak solution to (|3.ip with ^ G L^{VIt)- Then there exists a constant c, 
depending only on n, v, L, Li such that the following inequality holds whenever 
g(.To,to;2i?) C fl: 

(3.11) \u{xo,to)\ < cf {\u\ + s)dxdt + cll^\xo,tQ;2R) + cRs . 

jQ{xo.to-;R) 
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4. A NON-LOCAL CACCIOPPOLI'S INEQUALITY 

In [12 [12] we have developed more than one approach to the proof of the point- 
wise gradient estimates via non-linear potentials. Here we shall present one of these, 
taken form [22j . for the case p = 2, and for simplicity restricting to equations with 
no coefficients i.e. of the type 

(4.1) div a(i:'M) = ^. 

We believe that such method of proof is of independent technical interest since it 
potentially applies to all those problems with a lack of full differentiability, as it 
will be clear in a few lines. Moreover, we shall see that in the case (|4.ip estimate 
()2.7p holds component- wise; see ()4.1ip below. The assumptions considered for ()4.ip 
are of course 

(4.2) i^lAp < (a^(z)A, A) , \a,{z)\ < L , |a(0)|<L. 

which hold whenever z,A G R", where < ly < L. The presentation of this tech- 
nique is indeed one of the objectives of [22j . Aiming at the explanation of a general 
viewpoint, let us recall that for energy solutions u G VF^'^(ri) to homogeneous 
equations of the type 

(4.3) diva(i:»u) = 

the local boundedness of the gradient is achieved by first differentiating the equation 
(|4.3p . proving that Du G W^^^{fl), and then observing that v := D^^u solves the 
linear equation with measurable coefficients 

div{A{x)Dv) = A{x) a^{Du{x)) . 

At this stage the boundedness of D^v follows applying an iteration method, as for 
instance the one devised in the pioneering work of DeGiorgi [TT] . This is in turn 
based on the use of Caccioppoli's inequalities on level sets, that is, denoting 

(w — k)^ := max{w — fc, 0} , {w — k)^ :~ max{fc — w, 0} 

we have that inequalities of the type 

(4.4) / \D{D^u-k)+\^dx<^ [ \{D^u-k)+\'^dx 

JBr/^ ^ J Br, 2 

and similar variants, for instance involving {D^u — fc)_, hold whenever A: G K. In 
turn, the iteration of such inequalities yields the boundedness of D^u. In such 
an iteration, one controls the level sets of D^u via the higher order derivatives 
D{D^u — fc)-|- and Sobolev embedding theorem, building a geometric iteration in 
which, at every step, the gain is dictated by the Sobolev embedding exponent. 

Applying such a reasoning to the case (|4.ip seems to be difficult, as even in 
the simplest case (|l.ip it is in general false that Du G W^'^{n) when the right 
hand side /x is just a measure, or an L^-function. On the other hand, a result of 
[in] states that, although Calderon-Zygmund theory does not apply in the classical 
W^'^-sense, when considering the borderline case when ^ is a measure or lies in L^, 
it nevertheless holds provided the right functional setting is considered, i.e. using 
Fractional Sobolev spaces. Indeed, for SOLA to measure data problems as (|4.ip it 
holds that 

(4.5) Du G Wl'^'^n, R") for every e G (0, 1) , 

with related explicit a priori local estimates; see [THl Theorem 1.2] for precise state- 
ments. We here recall that, for a bounded open set A C M" and k € N, parameters 
a G (0, 1) and q G [I, oo), the fractional Sobolev space (^, M*^ ) consists of those 
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measurable mappings w: fl ^ MJ' such that the foUowing Gaghardo-type norm is 
finite: 

(4.6) =: + [w]a,q:A < oo . 

With such a notation (14.51) means that 



r„ 1 II \Du(x) ~ Du(y)\ , 



In' Jn' \x - y\ 

holds for every e € (0,1), and every subdomain Cl' d Q; the previous quantity is 
intuitively the L^-norm of the "(1 — e)-order derivative" of Du, roughly denotable 
by D^'^^Du. The inequality in (|4.7p let us think that Caccioppoli type inequality 
(|4.4p should be replaced by a fractional order version, and using the L^-norm, 
rather than the i^-one. Indeed we have the following theorem, that we again for 
simplicity state under the form of a priori estimate - i.e. assuming more regularity 
u £ W^''^{n) and fi £ L^(f2) (this can be again removed via an approximation 
scheme, and by considering suitable definitions of solutions). Needless to say, what 
it matters here is the precise form of the a priori estimate. 

Theorem 4.1 (Non-local Caccioppoli inequality). Let u G W^''^{n) be a weak 
solution to (|4.ip with fj. G L'^(fl), under the assumptions (|4.2p .' whenever ^ G 
{1, . . . ,n}, fc > 0, and whenever Bn C zs a ball with radius R, the inequality 

(4.8) [{\D,u\ - k)+U,BR,. " + £^i^KM , 

holds for every a < 1/2, where the constant c depends only on n, L, a . 



Comparing ()4.8p and (|4.4p , Theorem \J^.1\ tells us that for quasilinear equations 
CaccioppolVs inequalities are a robust tool that keeps holding at intermediate deriva- 
tives/integrability levels. We do think that the idea of using non-local Caccioppoli 
inequalities instead of the usual ones is interesting in itself as it leads to certain 
types of iterations which work without fully differentiating the equation; in turn, 
this could apply to all those problems with a lack of full differentiability. We indeed 
explicitly note here that a fractional Caccioppoli inequality has been indeed derived 
for notwithstanding the problems has integer order. The proof of the inequality is 
developed in [52] and has as a starting point some techniques introduced in [TH [TH] . 

The idea is now rather natural: inequality (|4.8p serves to start an iteration in 
which, at each stage we control the level set of D^u via the fractional derivative 
D^i^D^u) and the fractional version of Sobolev embedding theorem. We come up 
again with a geometric iteration whose step is in turn dictated by the fractional 
Sobolev embedding exponent. A point we want to emphasize, is that, as clearly 
inferrable from [52] , inequality (|4.8p contains all the information about the pointwise 
gradient estimate, no matter how small a is taken. As a matter of fact in the 
following we are not using explicitly the fact that m is a solution, but rather the 
fact that D^u satisfies (|4.8p . For this reason, we shall report the next result in an 
abstract way. Moreover, we think that the formulation below could be useful in 
different contexts. 

Theorem 4.2 (De Giorgi's fractional iteration). Let w G L^(f7) be a function with 
the property that there exist a G (0,1) and ci > 1, and a Radon measure fi, such 
that whenever Bn G n is a ball with radius R and k > 0, the inequality 

(4.9) [(i.,| - kuuusR,. ^ f X (i-i - + , 
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holds. Then the following estimate: 

(4.10) \w{xo)\<cJ \w\dx + cl[^^{xo,2R) 

Jb{xo,R} 

holds whenever B{xq,2R) C where the constant c depends on ci,n,a. 

The dependence of the constant c appearing in (|4.10p is not surprisingly as 
follows: 

lim c = oo and lim c = oo . 

Now we just have to conclude merging the last two theorems. Indeed, by Theorem 
14.11 we have assumption (|4.9p from Theorem 14.21 satisfied by w = D^u. In turn, 
applying Theorem l4.2l with such a choice of w we conclude with the desired pointwise 
gradient bound 

(4.11) \D^u(xo)\ < c f \D^u\dx + cl[^^{xo,2R). 

The last estimate clearly implies ()2.7p . being actually stronger since it holds for 
each single component of the gradient. 
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